Spin-Hall Effect in A Symmetric Quantum Wells by A Random Rashba Field 
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Changes dopant ion concentrations in the sides of a symmetric quantum well are known to create 
a random Rashba- type spin-orbit coupling. Here we demonstrate that, as a consequence, a finite size 
spin-Hall effect is also present. Our numerical algorithm estimates the result of the Kubo formula 
for the spin-Hall conductivity, by using a tight-binding approximation of the Hamiltonian in the 
framework of a time-dependent Green's function formalism, well suited for very large systems. 



PACS numbers: 71.70.Ej, 73.21.-b, 73.63.-b, 85.75.-d 

When first discussed almost fifty years ago^, the 
Rashba-type spin-orbit (SO) coupling was attributed to 
the inversion asymmetry in zinc-blende quantum wells 
grown along the [001] direction. In this picture, the SO 
interaction is linear in the electron momentum p and is 
written, in terms of the Pauli spin operators {cr^., CTj,, CTz} 
as: 



Hso = a{(Txky 



CTyh^ 



) 



(1) 



The coupling constant a is a constant proportional to the 
gradient of the electric potential across the well and thus 
tunable by an electric gate^. 

Searching for ways and means of manipulating the elec- 
tron spin in solid structures by electric fields, within the 
context of possible spintronics applications, has rekindled 
the interest in the Rashba model and numerous studies^i^ 
have been dedicated in the the past several years to un- 
derstanding all the implications of this interaction on 
electron transport in reduced dimensionality semiconduc- 
tor structures. Even though the original Rashba interac- 
tion is intrinsically linked to the quantum well asymme- 
try, a recent argument was made^ for the existence of 
a random SO-coupling that appears even in a perfectly 
symmetric quantum well on account of the changes in the 
dopant ion concentration on the sides of the well. This 
result is made possible by the existence of a local electric 
field, of random magnitude, perpendicular on the layer, 
at each point inside the well. Even though the spatial av- 
erage of the random Rashba field that ensues is zero, it 
has been demonstrated that it puts a certain imprint on 
various electronic properties, such as the spin-relaxation 
rate which acquires a minimum value. Moreover, in the 
presence of a magnetic field, the effects of the random 
Rashba fields lead to longer spin relaxation rates and a 
non-exponential spin-relaxation.^ 

These ideas suggest the existence of a finite spin- 
polarization that can be maintained also during spin 
transport. It is quite natural, therefore, to ask what hap- 



pens to the spin-Hall conductivity, in this situation. The 
present work investigates the effect of the random Rashba 
field in symmetric quantum wells on the spin-Hall con- 
ductivity, with the intention of appreciating its robust- 
ness against the natural variations of this type of SO cou- 
pling. Within a numerical algorithm based on the direct 
integration of the time-dependent modified Schrodinger 
equation, we establish that indeed, a minimum spin-Hall 
effect arises, dependent on the dopant ion concentration, 
that is quite resilient under the fluctuations of the ran- 
dom SO coupling. 

The physical model of our system consists of a two- 
dimensional (2D) quantum well sandwiched between two 
doped layers of different concentrations, separated by a 
distance zq. The dopant ions, assimilated with a (5-type 
perturbation, are assumed to have charge e and are local- 
ized by a 2D in-plane vector — {xi,yi). The total ion 
concentration in each layer is n(r) — J^j '^(i"~rj). In the 
simplest approximation the Rashba coupling constant is 
proportional with the z-component of the Coulomb elec- 
tric field generated by the impurities: 
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where e is the dielectric constant, and the summation is 
performed over all the dopant sites in both layers. Hence, 
a depends both on position r and on the dopant distri- 
bution. In general, we can write that Q;fl(r) = ancEzir), 
where is some phenomenological system-dependent 
constant parameter. An illustration of this analysis is 
presented in Fig [T] where values of aji{r) are showed for 
two different impurity concentrations, for a system size 
of 400 X 400nTO and for zq — bnm. 

Since, the position-dependent SO coupling constant 
aij(r) does no longer commute with the momentum op- 
erator, the single-particle Hamiltonian that describes the 
dynamics of an electron of momentum p and effective 
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FIG. 1: Spatial dependence of the Rashba spin-orbit inter- 
action for different values of Uu and Ud, the concentration 
of carriers on the up and down side of the well, respec- 
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and nd = 5 X 10 cm Right: 



tively. 

nd,u = 5 X lO^^cm"^. In both figures the system size is 
400 X 400nm and zo is fixed to 4nm. The energetic units 
are relative to to, the hopping integral. 



mass m* acquires a symmetric form 



1 



(3) 



where {A, B} represents the anticommutator of opera- 
tors A and B. The HamUtonian thus defined in Eq. ([3]) 
contains not just Hso^ the original Rashba Hamiltonian, 
Eq. I!]), written for aij(r), but also an additional part, 
Hrandom that reflects the non-commutativity of momen- 
tum and position operators, involved explicitly in aij(r). 



random 



dx 



dy 



(4) 



We note that Hrandom gives no contribution to other im- 
portant characteristics of the system, such as the spin- 
dependent force operator Fh ^ a|j(r), which is entirely 
generated by Hso^- 

Inspired by Eq. ([3]), we write the Hamiltonian of the 
total system of electrons directly within the tight bind- 
ing approximation, generating an expression suitable for 
numerical calculations: Thus, 



(5) 



<i,j>,a 



Here c]^ is the creation operator at site index i with spin 
a, Sx and Sy are unit vectors along the x and y direc- 
tions. An immediate consequence of the randomness of 
the Rashba spin-orbit strength is the site-dependence of 
the parameters that describe the model. First, an on-site 
energy, , that depends on the gradient components of 
the electric field appears, similar to the case of disor- 
der. The difference, however, is that in this case the 
disorder is spatially correlated. Then, the Rashba cou- 
pling Viiii) = Q;ij(i)/ao, written for a lattice constant ag, 
depends on the values of the electric fields at neighbor 



sites. The hopping coupling to has the usual expression 
to = l/2m*aQ. In the following estimates, we fix its value 
at to — 12meV which corresponds to a lattice constant 
ao — 2nm. (Throughout this work, h = 1, and free 
boundary conditions are used.) At the same time the 
charge concentration fixes the Fermi energy ep through 
the relation: n = 2mep/27rh^ = lO^^cm"^. A numer- 
ical estimate of the spin-Hall conductivity in a random 
environment, for large system sizes, can be obtained by 
using the Kubo formalism in the framework proposed by 
Tanaka^. This method can be, in principle, used for the 
calculation of any expectation value of any combinations 
of Green's functions and quantum operators, such as the 
density of states, or conductivity^!^. Besides, it is well 
suited for large systems where traditional methods such 
as direct diagonalization fails due to memory problems. 
Within this algorithm, for a given field distribution, the 
on-site energies, hopping probabilities and the Rashba 
field are computed first at each site. Then, the spin-Hall 
conductivity is calculated for a given configuration. 

The main steps involved in performing the numerical 
calculations are outlined below. The computation starts 
from solving the time-dependent modified Schrodinger 
equation, with a single-frequency source term: 



d |j,t) 



dt 



= i/|j,i) + |j)0(t)exp-''(^+''')* 



(6) 



where 77 is a finite small value and 9 is the step function. 

To determine the time evolution of the state ket | j; t), a 
direct numerical integration of the modified Schrodinger 
equation is performed, using the "leap-frog" algorithm^-. 
This is a second order, symmetrized differencing scheme, 
accurate up to (HAt)'^. In this approximation, Eq. ([6]) 
becomes: 

\j;t + At) ^ -2iAtH\j;t) + \j;t-At) (7) 
- 2iAt \j) e^'(^+"')*0(t) . 

The time step At is determined by At ~ P/\E\, where \E\ 
is the absolute value of the energy and /3 is a parameter 
whose value is less than 1 in order for the solution to be 
stable. 

An analytic solution of Eq. ^ , with the initial condi- 
tion |j, t = 0) = is written as 

|j,t) - -I f dt'e-'"^'-''^\j)e-'^''+'^^'' 
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\]) 
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where one recognizes the Green's function operator [E -\- 
irj — H)^^ as the prefactor in the final expression. For 
sufficiently large time, the Green's function operating on 
the ket |j) can be obtained with the relative accuracy 
5 = e^''-'" by inverting Eq. ([5]) as: 

G{E + iT?) \j) = lim |j, T) e*(^+*'')^. (9) 
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Consequently, the matrix element between any two states 
{i\ and |j) is obtained as: 



1 — >C>0 



(10) 



These results can be easily generalized to the case when 
a product including several Green's functions and oper- 
ators is involved by choosing a new initial state, such as 
= AG{E + irf) \i) in Eq. ([6]) and repeating the same 
procedure. In a similar fashion the method can be gen- 
eralized to the case when many different energy values 
are considered. In this case one solves Eq. ([B]) simul- 
taneously for a source term with multiple frequencies, 

The efhciency of the algorithm is increased if instead 
of a local orbital basis set it a randomized version of this 
basis is selected. This can be described described by a 
ket \4>) = '^n=i l*^) 6xp(—i (/)„), where \n) are the local- 
ized orbitals and (j)n are random numbers in the [0, 2tt] 
interval. Then, the average of a given operator A is cal- 
culated as 



\A\ 



J2in\A\n)+0i^ 



(11) 



within the statistical errors of 1/^/N. It was showedii 
that this choice of basis is the best one for reducing the 
numerical errors. 

As a first application of the present method we calcu- 
late the density of states (DOS), 



1 



3m [TrG (uj + irj)] 



(12) 



in a system size 200ao x 200ao for a constant Vr. The 
results, computed for three different values of Vr are pre- 
sented in Fig. O In the presence of a constant Rashba 
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FIG. 2: The density of states of a clean system with con- 
stant Rashba coupling Vr is calculated for a square lattice of 
size 200 X 200ao. The energy is expressed in hopping-integral 
units, to- 

field, the band structure suffers two important modifica- 
tions. First, the bandwidth increases proportional to Vr 



as a consequence of the renormalization of the hopping 
probability by the Rashba interaction. Second, the van- 
Hove singularity in the density of states that occurs at 
zero energy when Vr — splits when Vr, is turned on. 
This corresponds to the splitting and shifting of the orig- 
inal spin-degenerate band into the two chiral sub-bands, 
characteristic to the Rashba model . Analytically it can 
be shown^^ that the splitting is proportional with V^. 
This behavior can be observed in the inset of Figl2] where 
a large Vr = OMq was considered. 

The spin-Hall conductivity can be calculated with the 
Kubo formula 



X V, 



de 
2^ 
GA{e) - 



df{e) 



de 

fx GR[£)vy [GR{e 



{3l[GR{e)~GA{e)] 



(13) 



This expression considers contributions only from states 
at the Fermi surface, as indicated by the presence of the 
sharply peaked derivative of the distribution function. 
Such an approximation is justified by previous studies of 
the spin-Hall effect, which demonstrated that that states 
below the Fermi surface does not contribute to the spin- 
Hall conductivity and that spin-Hall conductivity comes 
entirely from quasiparticle states at the Fermi leveU^. 

In Eq. (fT3|) . the velocity operator is defined by the 
commutator: iVy = [y,H], whereas the spin current, 
along the x direction given in terms of the anticommu- 
tator between the velocity operator and the Pauli ma- 
trix cr^: j| — {a^jV^} /4^. Gr/a{() represents the re- 
tarded/advanced Green's function. 

For a constant Rashba field, the algorithm outlined 
above generates results, which after being averaged over 
200 different impurity distributions, are presented in Fig. 
131 For any value of Vr , the universally predicted value of 
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FIG. 3: Left: The spin Hall conductivity for a clean system 
of size 200 x 200ao, with a constant Vr. The results are 
obtained by averaging over 200 initial wave functions. Right: 
rj dependence of the spin-Hall conductance for a fixed energy 
E = 0.5t and a Rashba coupling Vn = Q.2t. The error bars 
are smaller than the symbols size. 

GsH — c/Stt is reached, but only for some energy inter- 
val inside the band. Usually the spin-hall conductivity 
is smaller than c/Stt. It decreases as the band edges are 
approached and vanishes beyond them. This behavior 
is preserved for system sizes up to 500ao x 500ao, in- 
dicating that the system remains in the ballistic regime 
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regardless of its size for as long as no disorder is included. 
In the right inset of Fig. Owe present the rj dependence 
of the spin-Hall conductivity. For values smaller than 
0.1 already the convergence to the approapriate value 
is reached. Increasing 77 above 0.1 deviations are al- 
ready consistent even if a larger integration time is used. 
Throughout of our calculations we have used i] = 0.1. 
Similar results for the spin-Hall conductance were ob- 
tained using a Landauer-Biittiker formalism^. 

Calculated values of the spin-Hall conductivity in the 
presence of a random Rashba field, averaged over 100 im- 
purity configurations, are showed in Fig. 0] for different 
dopant concentration of the sides of the well, u (up) and d 
(down) respectively. Here, the randomness of the electric 
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FIG. 4: The spin Hall conductivity -when both up-down dop- 
ing layers are present. For each Rashba field distribution 
the spin-Hall conductivity was averaged over 200 initial wave 
functions at each energy. The figure presents averages over 
100 random impurity configurations. 

field created by the dopants generate a Rashba coupling 
at each lattice site, even in the absence of an external 



Rashba field. A finite asH is still obtained. The parame- 
ters used for our analysis are experimentally reachable^. 
The solid and dashed-dotted lines in Fig. [4] correspond 
to the spatial distribution of the Rashba fields presented 
in Fig. Ill 

When a symmetric distribution of the impurity ions is 
considered on the sides of the well, a vanishingly small 
Rashba field ensues when averaged over the entire sam- 
ple, as shown in[TJ The corresponding spin-Hall conduc- 
tivity, however, which is determined by the configuration 
average of the spatial variation of the coupling constant 
a increases as expected with the impurity concentration 
and seems to be favored by perfectly symmetric distri- 
butions. When the numbers of impurities is different on 
the sides of the quantum well, the average Rashba field is 
finite, but the spin-Hall conductivity decreases. A possi- 
ble explanation of this result can be given in terms of the 
spatial correlation of the spin-orbit interaction over the 
spin precession length^ that are strong enough to assure 
the existence of a finite spin-Hall effect even in the case 
of a null average spin-orbit field. 

In conclusion we have studied the spin-Hall effect in 
a perfectly symmetric quantum well when fluctuations 
of the Rashba spin-orbit interaction field are considered 
as arising due to the impurities on the sides of the well. 
Even in the extreme limit, when an equal number of the 
impurities generate a Rashba field averaged to a very 
small value, a finite spin-hall conductivity seem to exist 
in the system in the range of 5 -;- 10% of the original 
universal value c/Stt. 
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